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SUMMARY 

Two-dimensional photoelastic stress analyses were conducted for epoxy 
resin models containing a hole defect under the conditions of Hertzian con- 
tact. Stress concentrations around the defect were determined as a function 
of several parameters. These were hole diameter, its vertical distance from 
the contact surface and the horizontal distance from the Hertzian contact 
area. Also determined was the effect of tangential traction (generated by a 
friction coefficient of 0.1) on the stress concentration. 

Sharp stress concentrations o®cur in the vicinity of both the left and 

the right side of the hole. The stress concentration becomes more distinct 

■ 

the larger the hole diameter and the smaller the distance between the hole 
and the contact surface. 
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The stress concentration is greatest when the disk imposing a normal 
load is located at the contact surface directly over the hole. The magni- 
tude and the location of stress concentration varies with the distance 
between the Hertzian contact area and the hole. Taking into account the 
stress amplitude, the area which can be involved in a process of rolling 
contact fatigue seems to be confined to a shallow region at both sides of 
the hole. The effect of tangential traction is comparatively small on the 
stress concentration around the hole. 

INTRODUCTION 

Although a great deal of effort has been expended in elucidating the 
mechanism of failures in concentrated rolling elements, the subject is still 
greatly debated, particularly with respect to what stresses are the critical 
ones causing failure. 

Most existing theories of contact fatigue damage are based on the 
classical Hertz theory and its subsurface related shear stresses. In the 

process of rolling contact fatigue, however, there are several important 

■ 

phenomena which cannot be explained satisfactorily by conventional con- 
siderations. Therefore, it has been pointed out for a long time that the 
mechanism must be understood in terms of several sources of stress concen- 
tration at or below the rolling contact surface. 

The numoer of papers analyzing the stress field, assuming the dis- 
continuity due to a defect such as a dent, an inclusion or a pre-existing 
crack in the contacting solid, is gradually increasing. There have been 
several studies of stress distribution analyses resulting from surface ir- 
regularities. Chiu and Liu presented a simple model for a furrow-shaped 
surface defect and analyzed the stress concentration around the defect 
(ref. 1). Leibensperger and Brittain assumed several kinds of asperity 
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ridge models and analyzed the micro-Hertzi an contact stress field in the 
vicinity of each asperity (ref. 2). 

Recently, several attempts were made to deal with the problem of crack 
initiation or propagation based on fracture mechanics. Keer and coworkers 
presented numerical results for a cracked half-space, surface- loaded by 
Hertzian contact stresses and introduced possible mechanisms for crack 
propagation (ref. 3). Hills and Ashelby tried stress analysis, paying 
attention to the development of subsurface cracks and discussed the stress 
involved in the propagation of the crack (ref. 4). 

There have been many experimental studies indicating that nonmetal lie 
oxide inclusions or voids contained in the material act as the predominant 
stress raisers for fatigue crack initiation. There is, however, little data 
concerning the stress distribution around the defects. Generally speaking, 
stress analyses of a contact problem while considering a discontinuity are 
highly complicated and, therefore, it is difficult to determine the stress 
distribution from take off theory. This study examines the simplest contact 

a 

problem assuming a circular hole-type defect in order to obtain the funda- 
mental data for a quantitative discussion of fatigue crack initiation. 

A conventional photoelastic technique was used, because it allows for 
the visualization of the stress concentration. In addition, the stress dis- 
tribution can easily be determined with an adequate model. The characteris- 
tics of stress distribution around the hole were analyzed with respect to 
several parameters. These were hole diameter, depth of the hole subsurface, 
position of Hertzian contact area and loading conditions. 

The results are discussed from the viewpoint of a mode of subsurface 
initiation of rolling contact fatigue. This is compared with the experi- 
mental results which have been reported concerning rolling contact phenomena. 
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STRESS ANALYSIS AND MODEL CONFIGURATION 

A conventional two-dimensional photoelastic technique was employed for 
stress analyses. A lens-system polariscope with polarizers and quarter wave 
plates 150-mm in diameter was set up across the epoxy resin model which was 
externally loaded by the disk. 

The size of the model was 110-mm in height, 175-mm in length and 6-mm 
in thickness, as indicated in figure 1. The photoelastic sensitivity of the 
epoxy resin used in this study is 1.00-mm/kg and its Young's modulus is 
280-kg/mm 2 . 

Models were prepared from rectangular shaped epoxy-resin plates 
(160-mm x 220-mm x 6-mm) available on the market for photoelastic stress 
analysis. The contact surface was finished longitudinally by emery paper 
about 0.5-um in peak-to-val ley roughness after milling. A circular hole 
simulating a defect of a real material was made at a specified location in 
the plate. The plate was drilled after being clamped between two dummy 
plates of the same material in order to produce a burr and crack-free hole 
edge. 

The hole diameter was chosen so that its ratio to the Hertzian contact 
width is as similar as possible to that frequently found in strength prob- 
lems of materials used in engineering applications. The lower limitation of 
the hole size is, however, determined inevitably by several practical re- 
strictions (i.e., tolerance in machining an exact shaped hole and photo- 
elastic resolution). Therefore, the hole diameter was selected to be within 
tne range of 1 to 3-mm. 

Loading 

An external load was applied to the model by a smooth surface finished 
disk made of acrylic acid resin (300-mm in diameter and 8-mm in thickness) 
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using the loading system snown in figure 1. Normal and tangential load can 
be applied by lever arrangements. In most experiments a normal load of 
77— kg was applied, taking into account the mechanical strength of contacting 
materials and the extent of the isochromatic fringe pattern in the model. 

The Hertzian contact width was found to be approximately 7-mm in these 
cases. For some experiments, a disk of polycarbonate resin (5-mm in thick- 
ness) was used. Except for a slight difference in Hertzian contact width, 
no difference can be observed in the stress distribution characteristics 
obtained with the two disks. 

Analysis 

Prior to each experiment the model was annealed and cured in a 
temperature-control led furnace in order to eliminate the residual stresses 
brought about during the model-making process and by the so-called time-edge 
effect. 

In photoelastic analyses, two kinds of information can usually be 
obtained: (1) the isochromatics, or loci of points at which the value of 

maximum principal shear stress is constant, and (2) the isoclinics, or loci 
of points which exhibit the same inclination angle of the principal 
stresses. It was difficult, however, to obtain fine and resolvable iso- 
clinics corresponding to the highly ordered isochromatic fringes which can 
be observed even with a simple polariscopic system and a reasonable model 
size. Therefore, in most analyses the distribution of maximum shear stress, 
that is, a half difference in principal stresses, was determined exclusively 
by paying attention to the isochromatic pattern around the hold. At a free 
boundary of a peripheral surface of the hole along which no external forces 
are applied, one of the principal stresses normal to the surface is zero. 
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Therefore, the remaining principal stress parallel to the boundary surface 
can be determined uniquely by the isochromatic fringe order. 

Analytical results obtained in such a manner could, therefore, offer 
information concerning the mechanism of rolling contact fatigue. 

Isochromatic photographs were taken in a dark background (i.e., in a 
circular polarized light polariscope) . The nomenclature employed in this 
study is shown in figure 2. 

RESULTS 

Representative isochromatic patterns beneath the disk are shown in fig- 
ure 3. These were obtained with a hole and without a hole. The stress- 
raising effect of the hole appears distinctly, since fringes having the num- 
ber of higher order can be recognized at the area close to the hole edge 
(fig. 3(b)). The stress concentration around the hole is indicated in fig- 
ure 4. A circumferential principal stress located at a point on the hole 
edge can be determined from the number of ordered isochromatic fringes which 
intersect with the hole edge at the same location (fig. 4(a)). The stress 
distribution determined by isochromatic fringes around the hole is indicated 
with a solid line in figure 4(b). The dotted line represents the stress 
distribution estimated from a theoretical analysis by Tamate for the case of 
a contact between a rigid stamp having a flat surface of finite width and a 
half-spaced elastic solid (ref. 5). 

Certainly there are marked differences between the two analyses. 
Pnotoelastic fringe patterns do not, however, offer a substantially differ- 
ent pattern between the two, especially around the hole (fig. 5). There- 
fore, it is reasonable to compare the results directly with each other. 

Both of the results in figure 4(b) show a similar pattern in the stress 
distribution. A large amount of stress concentration is observed at both 
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the left and the right side of the hole, the location of which is denoted at 

tne -90 and +90 degree points, respectively, as was defined in figure 2. 

The magnitude of the stress located at these regions reaches 1.5 P max as 

determined by tne photoelastic analysis, where P max means tne maximum 

Hertzian pressure obtained under the condition of no hole. The maximum 

shear stress t , , that is, a half-circumferential stress in this study, 
max 

reaches 0.75 P m3 , which exceeds twice the maximum shear stress 
max 

0.30 P found to occur beneath the contact surface having no hole, 
max 

The stress at the same location obtained from a theoretical estimation 
is, however, larger. The difference may be caused either by differences in 
loading conditions or configurations. 

A theoretical estimation shows a slight tensile stress appeared near 
the bottom of the hole, while photoelastic analysis indicates all the stress 
is compressive around the hole. Even if it is small in magnitude, a tensile 
stress may act detrimentally in the fracture process of practical components 
under conditions of Hertzian contact. Additional analyses may be necessary 
in order to determine the stress distribution in more detail at the area 
where a tensile stress is predicted to occur. Photoelastic resolution in 
this study was not satisfactory at the upper and lower areas of the hole. 

In this paper, analytical conditions were so taken as to cause no ten- 
sile stress at tne peripheral surface. Since the maximum shear stress has 
been regarded as a primary stress causing various characteristic phenomena 
concerning the rolling contact failure, much attention was paid exclusively 
to the stress distribution at the left and the right side of the hole where 
a relatively large amount of shear stress concentration occurs. 

The analytical results presented herein are represented by the term 
T max ^ or max ^ mum shear stress as a matter of convenience. In addition. 
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attention was paid to the stress distribution at the right half of defects, 
since defects are hole-shaped and symmetrical. 

The stress concentration around the hole obtained by varying the hole 
diameter are represented in figure 6. Figure 6(b) shows their stress dis- 
tributions and figure 6(a) indicates both variations of the greatest maximum 
★ ★ 

shear stress r and its location e on the hole edge with hole 
max 3 

* * 

diameter. The stress x „ increases, but its location o scarce- 

L max 

★ 

ly changes with the magnitude of hole diameter. As expected, x max 

finally reaches zero with an extrapolation of the decrease of hole diameter, 

which is plotted in the figure by a dotted line. 

Figure 7 presents the relationship between the stress concentration and 

the distance of the hole from the contact surface. Stress distributions are 

★ ★ 

indicated in figure 7(b) and variations of Tmax and o are 

★ 

plotted against the distance. The stress x max decreases uniformly 

★ 

with an increase in the depth of hole and its location e changes 
slightly within a range of 10 degrees by varying the non-dimensional depth 

a 

over the range of 0.50 to 1.5 d/b. 

It is well known that tangential traction located at the contact sur- 
face has a significant influence on rolling contact fatigue (ref. 6). An 
example of the effect of tangential load on the stress distribution is shown 
in figure 8. The radial ordinate is graduated in numbers of fringe order so 
as to magnify the variation of the effect. The magnitude of maximum shear 
stress, x m , v , is one-twelfth of the number of the fringe order N, that 

ilia A 

2 

is, i max = 1/12 N (kg/mm ). Therefore, the effect of tangential load is 
relatively small, although the coefficient of tangential traction, 0.1, is 
larger than that usually predicted in practical rolling contact applica- 

t 
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tions. The effect of this load on the stress concentration around a hole- 


type defect is, therefore, concluded to be negligibly small. 

Components in rolling contact applications usually experience repeti- 
tive loading. Therefore, in order to discuss the mechanisms of rolling con- 
tact fatigue, it is necessary to obtain the stress variation over which an 
element in a component is exposed during the course of rolling contact. 

The stress distribution behavior changes with the variation of horizon- 
tal distance between the hole and the Hertzian contact area as is illus- 
trated in figure 9. The photograph of figure 9 shows a representative iso- 
chromatic pattern obtained under the contact condition in which the Hertzian 
contact area deviates from the location directly over the hole. Locations 
where marked stress concentrations are observed shift along the hole edge, 
depending on the horizontal distance. The magnitude of the largest maximum 
shear stress decreases and its location deviates from the ±90 degree point 
with increasing distance. The variations on the right half of the hole edge 
plotted against the distance under several conditions are represented in 
figure 10. Although it has already been pointed out in the descriptions of 
the individual effects, hole size, location of the hole from the surface and 
tangential Toad do not change dramatically in the characteristics of their 
behavior by varying the position of Hertzian contact area. 

These relationships do not give exact stress amplitudes in certain 
fixed directions. They are, however, useful for discussing where an incipi- 
ent micro-crack initiates, because an area which experiences a severe stress 
concentration repeatedly can be considered as one of principal sources for 
crack formation during the rolling contact process. 
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DISCUSSION 


Recently, a number of reports have indicated that a large number of all 
failures in rolling contact elements is associated with severe surface de- 
fects or stress raisers. 

However, Lundberg and Palmgren's type of life prediction formulation 
has been frequently quoted and applied to practical rolling contact prob- 
lems. This prediction was developed assuming that most failures start in 
the subsurface zone of high shear stress below the rolling contact surface. 
This situation seems to occur since there is ample evidence which indicates 
that several metallurgical phenomena, deformations or micro-structural 
alterations, for example, take place in the subsurface high-stress zone 
(refs. 7 to 9). Furthermore, micro-cracks can be frequently observed 
propagating near or along this area, the metal lurgical character of which is 
different from the matrix. It is still debated whether or not the deformed 
area and stress-altered micro-structures are: (1) a primary source for 

initiation of micro-cracks leading to the final failure, ( 2 ) these areas 

» 

facilitate the propagation of cracks that are initiated elsewhere, or (3) do 
not take a significant part in the rolling contact failure. 

It is well known, however, based on the experimental results, that 
fatigue life is shortened by pre-existing defects such as nonmetallic in- 
clusions, voids or other metallurgical defects. 

In this study, a simple hole-type defect was assumed existing in the 
matrix and a large amount of stress concentration could be observed in the 
vicinity of tne hole as expected. The stress concentration could be ana- 
lyzed accurately. The magnitude of the stress concentration was observed to 
increase with increasing hole diameter. Conversely, the smaller the hole 
size, tne less significant is the influence of the defect. Figure 6(a) 
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suggests that there is a critical size below which the magnitude of stress 
concentration is smaller than the ordinary maximum shear stress 
(0.301 P ) occurring under the condition of Hertzian contact and, there- 
fore, the influence of the defect becomes negligible, when the location of 
the hole is shifted up close to the surface, the magnitude of the stress 
concentration increases uniformly and accordingly a more detrimental in- 
fluence might be produced on the rolling contact failure. 

If the hole diameter exceeds the size examined herein, a more intensi- 
fied stress concentration might be obtained. In such a case, however, the 
ratio of the hole diameter to the Hertzian contact width departs from a 
reasonable value encountered in practical problems. Further, the tensile 
stress is under such circumstances predicted to occur near the upper side of 
the hole. Therefore, the characteristics of such a problem are substan- 
tially different from those of the present study. 

The results obtained in this study are similar to the experimental re- 
sults discussed by Soda for the effect of inclusions existing in materials 
(ref. 10). These results also indicate that even when the defect is located 
near the classical high-shear stress zone (0.786b), the stress concentration 
is not emphasized in the sense that there is an interaction between the 
stress field due to the defect and the stress field associated with the 
ordinary Hertzian contact. Therefore, if the defects located near the sub- 
surface nigh-stress zone have a harmful influence on the rolling contact 
fatigue, this effect is better understood with the mechanism increasing the 
probability for the occurrence of failure due to the superposition of the 
defect and conventional high-shear stress zone, and not by emphasizing the 
stress concentration due to the interaction of these two effects. 
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On the other hand, many investigators have reported cracks propagating 
parallel to the rolling contact surface in both practical applications and 
bench tests. In the case of a material containing inclusions, incipient 
micro-cracks have been frequently observed to the right or left side of 
their peripheral surfaces. 

In a more exact discussion of the mechanism of rolling contact fatigue, 
it is necessary to formulate the fracture criterion applicable to the roll- 
ing contact fatigue and to determine all the stress components contributing 
to the crack formation in the fatigue process. 

The results of this study, however, show that even when the Hertzian 
contact area deviates from the hole, the stress concentration takes place at 
the restricted area near the left and right side of the hole edge. Since 
the maximum shear stress, x max decreases with the deviation, the area 
where possible damage seems to accumulate causing a micro-crack should 
further contract toward the ±90-degree position. Therefore, the results 
correspond well with rolling contact fatigue cracks observed in the vicinity 

a 

of inclusions. 

The contours of stress distributions around the hole are very similar 
morphologically to the so called "butterfly" which can be observed near an 
inclusion in a bearing material undergoing cyclic stresses. The analyses 
should be extended to include a discussion of the mechanism of formation of 
the "butterfly." Since a crack initiated in a material may change the 
stress distribution, alternate analyses must be conducted to explain crack 
propagation. 

The magnitude of the tangential traction, in this study, is larger than 
that encountered in practical systems. However, the effect of the tangen- 
tial traction on the stress concentration is relatively small. This phe- 
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nomenon predicts that in the case of incipient microcracks initiating 
beneath the surface, the tangential traction does not have a significant 
role. This tendency should be compared with the fact that for the surface- 
originated mode of rolling contact fatigue there is a marked influence of 
tangential traction (ref. 6). 

CONCLUSIONS 

The characteristics of the stress distribution around a hole-type 
defect were described under the condition of Hertzian contact with respect 
to several parameters; that is, loading, hole diameter and the location of 
tne defect (hole). The following conclusions are drawn from the results 
obtained: 

1. A high stress concentration, which has a high stress gradient, 
exceeding the conventional maximum shear stress, can be observed in the 
vicinity of the defect (hole). Analytical results indicate a similar tend- 
ency to that observed in rolling contact fatigue studies for materials con- 
taining inclusions. 

2. Taking into account a stress amplitude which is brought about by 
varying the distance between the position of the defect (hole) and contact 
area, the greatest deterioration may be concentrated at a restricted area 
near the left and right hand sides of the peripheral surface of the defect. 

3. Tangential traction has only a slight influence on the stress con- 
centration around the defect (hole). 

Professor T. Hirata (Tokyo University of Agriculture and Technology) 
allowed us to use the polariscope. 
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(a) ORDINARY ISOCHROMATIC PATTERN BELOW DISK (NO HOLE). 



(b) ISOCHROMATIC PATTERN OF THE MODEL CONTAINING A HOLE - 
2a = 1.0 mm (0. 29 b), d = 3. 0 mm (0. 88 b). 


Figure 3. - Representative isochromatic patterns below disk - normal 
load = 61kg, P^ = 1.90kg/mm^, 2b = 6.8 mm. Number in the 
figure denotes isochromatic fringe order. 
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CIRCUMFERENTIAL STRESS, 
a 0 /P max 



(b) 


Figure 4. - Isochromatic pattern and circumfer- 
ential stress distribution. 

(a) Isochromatic fringes around the hole and fringe 
order; 2a = 1. 0 mm, d = 3. 0 mm, 

normal load = 77 kg. P max = Z 33 kg/mm 2 . 

(b) Circumferential stress distributions around the hole-, 

Solid line: Result by photoelasticity (in the case of 
contact against disk) b/d = 1. 3, a/d = 0. 2. 

Dotted line: Estimation by theory (in the case of 
contact against flat-ended rigid stamp) b/d = 1. 0, 
a/d - 0. Z 
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NORMAL LOAD 
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■ 67 kg, 


2a = 2. 0 mm, 2b = 6. 0 mm, 
d = 3. 0 mm 



(a) (b) 


Figure 5. - Isochromatic patterns below disk and rectangular stamp. 

(a) Isochromatic pattern below disk; b/d = 1. 3, a/d = 0. 3. 

(b) Isochromatic pattern below flat-ended rigid stamp; b/d = 1. 0, a/d = 0. 3. 
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Figure 6. - Effect of hole diameter on stress concentra- 
tion ; d = 3. 0 mm (d/b = 0. 79), p= 0, E; = 0. 

(a) Variations of the greatest maximum shear stress and 
its location with hole diameter. 

(b) Maximum shear stress distributions around holes. 
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Figure 7. - Effect of distance between hole and contact 
surface on stress concentration; 2a = 1.0 mm 
(a/b = 0. 13). H - 0. E; - 0. 

(a) Variations of the greatest maximum shear stress and 

its location with distance between hole and contact surface. 

(b) Maximum shear stress distributions around holes. 
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FRINGE ORDER 



(c) 

Figure 8. - Example of effect of tangential force on stress concentration; 2a = 2 mm 
(a/b - 0. 26), d = 3 mm td/b = 0.79). 

(a) Isochromatic pattern; p= 0.1. 

(b) Isochromatic pattern; p= 0. 

(c) Stress distributions around the hole. 

• In the case of (a) 

O In the case of (b) 
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MAXIMUM SHEAR STRESS 
T max^max 



(a) 



(c) 



Figure 9. - Effect of distance between hole and Hertzian contact area on stress con- 
centration; 2a = 2. 0 mm (a/b ■ 0. 26), d « 3. 0 mm (dfb « 0. 79), p* 0. 

(a) E/b - -1.5; (b) E/b - -1. 3 ; (c) E/b - -0.5 ; (d) E/b * 0.2; (e) E/b ■ 1.5. 
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o 2a * 1. 0 mm (a/b = 0. 13), d = 1. 9 mm 
(d/b = 0.5), p=0 

• 2a = 1. 0 mm (a/b = 0. 13), d = 1. 9 mm 
(d/b = 0. 5), (j= 0. 1 
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Figure 10. - Variations of the greatest maximum shear 
stress and its location with distance between hole and 
Hertzian contact area under various conditions. 
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